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ABSTRACT 

Adaptive arrays are highly versatile sensors employed in modem wireless communication 
systems to combat interference and multi-path fading and thereby increase system capacity. In 
adaptive processing, weights are attached to the incoming signal at each element to produce nulls 
in the directions of interferers. However, mutual coupling is normally ignored in such processing. 
Instead, the principle of pattern multiplication is used, where the assumption is that the radiation 
pattern of an array is the individual radiation pattern of the elements multiplied by an array 
factor. This assumption ignores mutual coupling and the error can be significant under certain 
conditions. This work begins by setting up the scenario of an adaptive array in a mobile 
communications scenario in a scattering environment. Following that, we introduce the theories 
of mutual coupling and beam-forming. Expressions for the optimum solutions and the signal-to- 
interference-plus-noise ratio are then derived based on the preceding discussions. Subsequently, 
we show the effects of mutual coupling and the conditions under which the effect is significant. 
In addition, we will examine various parameters involved in our scenario and how they affect the 
performance of the adaptive array. 
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I. INTRODUCTION 



Array antennas are well known for their versatility in many areas, not least of 
which in point-to-point communication. By a process of beam-forming, an array antenna 
can steer its main beam in the direction of its intended transmitter or receiver while 
simultaneously pointing its nulls to cancel interference sources. Furthermore, adaptive 
algorithms have been introduced [Ref. 1 : p586 - 591] to infer and respond to the 
directions of these (desired and interference) sources, the positions of which change with 
time. There is currently a lot of interest in using adaptive antennas in mobile 
communications to increase systems capacity or reduce interference. 

However, the effect of mutual coupling has received little attention in the adaptive 
array processing community. It is often assumed that the beam pattern of an array is the 
individual pattern of an element multiplied by an array factor, in which mutual coupling 
is ignored. In this work, we will show that this error is significant under certain 
conditions and examine the effect mutual coupling has on the performance of array 
antennas. 

We begin in Chapter II by describing the scenario of mobile communications on 
which this work is based, followed by theoretical discussions on mutual coupling and 
beam-forming. Finally, we shall derive useful expressions for optimum adaptive solutions 
and for the signal-to-interference-plus-noise ratio (SINR), which will serve as the 
yardstick for measuring the performance of an array antenna. 

In Chapter m, we will show the effect of mutual coupling and the conditions 
under which this effect is significant. In addition, we will examine various parameters 
involved in our scenario and how they affect the performance of the adaptive array. 
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II. THEORY 



In this work we will look at the effect of element mutual coupling on the signal- 
to-interference-plus-noise ratio (SINR) of an adaptive array whose design is based on the 
minimum mean square error algorithm. The elements of the array are assumed to be 
dipoles and the analysis will be carried out in terms of the array mutual impedance. A 
variable number of sources of interference as well as variable inter-element spacing will 
be considered in the study. 

A. DESCRIPTION OF THE PROBLEM 

We consider the case of K+l transmitters (one desired and K interference) all 
radiating narrowband signals to a receiver at the same carrier frequency. The scenario is 
illustrated in Fig 2.1. The transmitters (also referred to as sources) are all embedded in a 
scattering environment. The scatterers, which could represent reflecting obstacles 
surrounding the source, are assumed to be distributed uniformly around each source over 
a circular region of radius r^ We assume single bounce scattering of waves, meaning that 
waves reach the receiver via single scattering off the obstacles. We further assume a 2D 
situation where waves travel only in the azimuthal (horizontal) plane. Because of the 
presence of scatterers, waves will not arrive just from a single direction for each source 
but from an angular region determined by rk and the distance Rk between the k th source 
and the receiver. The receiving antenna is assumed to be an N-element linear array with 
an inter-element spacing d. The elements are assumed to be center-fed half wave dipoles 
oriented along the z axis, so the pattern of an isolated element in the azimuth plane is 
constant. The mean directions of the desired source (k=0) and the k th interferer (k=l to K) 



are denoted by (j ) 0 and ^respectively. Let (j) k be the angle from any scatterer around the k th 
source to the array and - <p k , as shown. The angle of arrival £% is then measured 

from the mean direction <j> k for each source. 




Because of the multi-path phenomenon, the amplitude, phase and angle of arrival 
of waves from each source will be random variables. The interferes are assumed to be 
arbitrarily located about the receiver. The purpose of the array is two-fold: 

i) the array weights are chosen to produce nulls in the direction of the 
interferes, and thereby reduce their influence on the received signal, and 

ii) if N > K+l, to provide diversity gain against multi-path fading. 
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The weights of the array are continually changed depending on the location of the 
interferes. Such an array is termed as adaptive array or smart antenna. The situation 
considered above may well pertain to the case of wireless network in a tactical battlefield 
or to the case of mobile cellular radio where K+l co-channel mobiles send signals to a 
base-station receiver. We have assumed thus far that the mean directions of the desired 
source and the interferes are known a priori. This is seldom the case in reality. In a 
practical mobile situation, the sources move around in space and their directions are not 
known precisely. In such a case, the weights are determined by means of either the Least 
Mean Square (LMS) or Recursive Least Square (RLS) algorithm. However the present 
case is still of interest in applications such as wireless local loops where wireless service 
is provided to fixed subscribers and the directions of subscribers available from a 
previously compiled data base. 

Since the calculation of optimum weights requires the covariance matrix of the 
received signal vector, we would like to know the probability density function (pdf) of the 
angle of arrival. For a uniform density of scatterers about the source, the pdf of the angle 
of arrival with respect to the mean direction is given as [Ref. 2 : p24-25] : 

p(otk) = :1 - L cosa k jl-b k 2 sin 2 a k |od<£* (1) 

K 

where bk = Rk / r k (2) 

and ^k is the maximum angle of arrival from the k 111 scattering region, given by 

bkSin^k = 1 

For subsequent analysis, it is simpler to express the pdf in terms of a new angle 7k 
defined by siny k = bksinctk ; then 
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cos 7k ctyk = b k cos a k da k 



(3) 



Using this transformation of variable [Ref. 3 : pi 22], the pdf of Yk is: 



P(7k) = P(OCk) 



da k 



dYk 



2 2 
-cos Y k 
n 



K ^ 

-Yk- — 

2 2 



(4) 



This will be applied in Section ID when the expressions for optimum weights and SINR 
are derived. 
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The mean value of a k is a k = 0 and its variance oc£ for b k »l is approximately 



given as [Ref. 4] : 



« = C i- 



2 3 



k < 



(5) 



\2 

1 

V 2b k y 

v K y 



2 j 



r 

The root mean square (rms) value of the angular spread is then . For example, with r k 



= 100 m and R k = 1 km, ^ k = 5.7° and a. =2.85°. 



B. MUTUAL COUPLING 

In an array of electrically identical antenna elements, it is common to assume the 
total radiation pattern to be the product of the array factor and the radiation pattern of the 
individual elements, i.e., total radiation pattern = (array factor) x (element radiation 
pattern). Implicit in this formulation is the assumption that the individual radiation 
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patterns remain the same even when placed in proximity of each other when mutual 
electromagnetic interaction exists. However, such an assumption is erroneous in view of 
element mutual coupling. When an antenna array is designed to provide diversity gain, 
mutual coupling will disturb the spatial correlation of signals received by the array. In this 
thesis, we will examine the effect of mutual coupling on the performance of the array. 

The input impedance of an isolated antenna is simply the ratio of the voltage 
applied (at its terminals) to the resulting terminal current. This is known as the self- 
impedance and is, in general, a complex value as both voltage and current are of different 
phases. When two antennas are located in close proximity, the current distribution (and 
hence the radiation pattern) on one is affected by the field radiated from the other. This is 
known as mutual coupling. To account for the mutual coupling effect, extra impedance 
terms called mutual impedances are needed. 

In general, the solution to the mutual impedance for dipoles can be accomplished 
by way of coupled integral equations. Apart from its dependence on physical dimensions 
and separation distance between the elements, the mutual impedance also changes with 
the number of antenna elements. For an array of N elements, the impedance matrix can be 
written as Znm [N ^ denoting the mutual impedance on the m th element induced by the n 1 * 1 
element in an N-element array. The details of mutual impedance calculation are not 
within the scope of this paper; we are interested in the effects it has on array antennas. We 
will use results from simulation programs or from analytical expressions. As mentioned 
earlier, for the sake of simplicity, we shall only concern ourselves with half wavelength 
dipoles of small thickness. Even though Zn, n [N] varies with N, we will assume Znm tN1 = 
Zmn 121 ; actual computation of Zm n [4] and Zm,/ 81 show them to be very close to Znm t2] . (The 
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superscript [2J is subsequently dropped.) Z™, the mutual impedance between two side-by- 
side dipoles, each of length L and separated by a distance d, is given by [Ref. 5: p265] : 

Z 12 = 30{ 2Ej (- jk 0 d) - Ej [-jk 0 (Vd 2 +L 2 + L)] - Ej [- jk 0 (Vd 2 + L 2 - L)] } (6) 

where k 0 = 27C/A, is the wavenumber in free space, and 
E,(±jy) = Ci(y)± jSj(y) 

Q(y) = cosine integral 
Sj(y) = sine integral 

It is obvious that the inter-element spacing d has an influence on Z™,; as the 
distance apart increases, one can expect the contribution from other array elements to be 
reduced and hence Z™ (m ^n) diminishes. The resultant impedance matrix Z is then 
almost diagonal and mutual coupling can be ignored. 

C. BEAM-FORMING ALGORITHM 

Consider an arrangement of N elements spaced d apart (Fig 2.2). An 
electromagnetic wave arriving from bearing (|) induces an open circuited voltage V n oc on 
the n th element. The open circuited voltage depends on the incident field and the 
orientation of the dipoles with respect to the incoming wave. The coupling between the 
elements is governed by the impedance matrix Z. It is assumed that the adaptive 
processor presents load impedance (ZOn (n=l, 2, ... N) at each element. The terminal 
voltage and current at the n th port are denoted by V n and I n respectively. It is clear that V= 
-Z L I, where V is a vector of voltages presented to the adaptive processor and Zl is the 
diagonal matrix = diag(Z L i , Z L 2 , • . . , Z L n). 
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Fig 2.2 Array of N elements 

The terminal voltage is related to the open circuited voltage by means of the 
impedance matrix as : 

V, = Z n Ii + Z 12 I 2 + Z 13 I 3 + - . + Z 1 N I N + V, oc 
V 2 = Z 2 ,Ii + Z 22 I 2 + Z 23 I 3 + ... + Z 2 nIn + V 2 ° c 
V 3 = Z 3 ,I, + Z 32 I 2 + Z 33 I 3 + ... + Z 3 nIn + V 3 0C 



Vn = ZniIi + Zn 2 I 2 + Zn 3 I 3 + . . . + ZnnIn + Vn° c 
or equivalently, 

V = ZI + V oc (7) 

For a dipole antenna the open circuited voltage can be obtained from the effective 
length of the antenna by 
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(8) 



V n oc = h.E“ 



inc 



where 



h = vector effective length of n th element 



Ep 10 = incident electric field on n th element. 

For dipoles of length L, oriented along the z-axis and having a sinusoidal 
distribution I(z) = Io sin { k 0 (L/2-lzl) } , the effective length is 

X f(0) 



h = e 



n . ,k 0 L 
s m (— ) 



(9) 



where 0 is the usual polar coordinate in the spherical coordinate system and 



f(0) 



k n L ^ k^L 
cos( — - — cos 0) - cos( — ^ — ) 

sin0 



is the radiation pattern of the dipole. 

The electric field incident on the n th dipole is of the form, 
ginc_ gg^ e jk 0 d(n-l)cos(t> 

In the present case, 0 = n / 2 and the radiation pattern of the dipole reduces to 



f(0) = (l - cos ^L) 



, 1 — COS(— ^— ) 

Therefore, V n oc = -E 0 2 e jk 0 d(n-l)cos<|> = Voe jk 0 d(n-l)cos<l) (10 ) 

s,n< — ) 

where V 0 = — - tan( — - ) 

K 4 

Using (10) in (7) and substituting V = -Zl I , we get 

V = Z C V 0C (11) 
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where 


Z c = Zl Z 0 1 


(12) 


and 


Zo = (Z L + Z) 


(13) 



It is this voltage vector V and not V oc that is available for array processing. It is 
seen that the effect of mutual coupling is to couple the open circuited voltages induced at 
each element. Thus mutual coupling will affect the spatial correlation between the 
element signals and hence affect the diversity gain of the antenna. 

In the presence of thermal noise, the total voltage input to the adaptive processor 
is 

V, = Z c V oc + N(t) (14) 

where N(t) = column vector of thermal noise. 

With this input signal vector V t , the adaptive processor attaches a weight w n * (* is 
the complex conjugate operator) to each (V t ) n as shown in Fig 2.3, and produces a 
combined output S r . 



1 


(V t ), 

r i 


(V,) 2 

r 




(V,) N 

r 


1 


p— «-wi C 

f L 1 


P — 4-w 2 * 

- ► 


, « 1 


p— ◄^VN* 




W V 








1 


r 

Sr 



Fig 2.3 Adaptive Processor 



In mathematical terms, the array response is 

S r =wi Vi + W2* V2 +.... + wn* Vn 

= w H V t (15) 
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where H denotes the Hermitian operator (transpose and conjugate). 

By varying the weighting vector w, up to N-l nulls can be steered in any direction. 
For example, in the absence of mutual coupling and scattering, and using a 4-element 
array, if the nulls are desired at <J>x , ((b and (fo we choose the weights according to 



(V, 


v, 2 




w, 




r-n 


V 2 




Vf 




= 


-i 


V, 


V 3 2 


V , 3 


W l 




-i 


l 3 




3 ) 


w 4 




k j 








V w - J 







where V k = e jk ° d(n 1)cos<t>k , k = 1,2,3. For <)>, = 45°, <t> 2 = 90°, 4> 3 = 135° and d = X/2, 
W 7 Xi\ w d 

we get — = 0.2114, — = -0.2114 and — = 1. The array response to a plane wave 
w, w, w, 

arriving from angle (j) is plotted in Fig 2.4. 



4 elements, each 0.5 X apart. 




Fig 2.4 Beam pattern of a 4 element array, intended nulls at (J) = 45°, 90° and 135° 
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1 . 



Criteria for Optimum Weight 



It was shown in the preceding paragraphs that, by using an appropriate weighting 
vector, one is able to direct the nulls towards the direction of interference sources, thereby 
cancelling them out. How does one choose the weights in the presence of noise and 
scatterers? 

The choice of the weighting function is based on the statistics of the total signal 
received (desired + interference). Basically, the objective is to optimize the beam towards 
the desired signal according to a prescribed criterion, of which there are several. It has 
been shown [Ref. 5: p41] that the various criteria (including the Minimum Mean Square 
Error (MMSE) and the Maximum SINR criteria) converge to the same Wiener solution, 
differing only by a trivial scalar factor. Since the solutions are all the same, we shall use 
the MMSE criterion. In the MMSE algorithm, the weights are chosen to minimize the 
error, defined as: 



where < > is the expectation operator. The weights that minimize the error are shown to 
be of the form [Ref. 6: p39] : 



error = < | A r (t) - S r (t) | 2 > 



(17) 



w opt = O S 



(18) 



where 



= covariance matrix = <V t V t H > 



(19) 



S = <A r H V t > 



( 20 ) 



A r = reference signal = aoexp[j(cot + (3 0 )] 
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2. Choice of inter-element spacing d 

Before proceeding to the calculation of the weight vectors, it would be instructive 
to consider the effect the separation distance between array elements on the input signal. 
Suppose there were no mutual coupling, no thermal noise (N(t)=0) and all loads are 

2 

identical (hence Z c = z c U , where U is an NxN identity matrix and z c = ). Then 

Z L +z ll 



using (10) - (14), the input signal to the array processor is : 
(V,) n = z c V 0 e-> k ° d(n - 1)cos * 



( 21 ) 



Suppose the weights are all equal, i.e. w n = w for all n, then 

s,«>) = i>„*v„ 



1 



= w*z c V 0 ^ k » d(D - 1)cos ^ 
1 



1 _ e j- Nk o d COS <t> 

w * z r V n 

j _ e J k o d cos f 



The maxima are obtained at cos <j) = ± n ( 




( 22 ) 



(n=0,l,2...) where the denominator 



becomes zero. In general, there is the desired major lobe and other major lobes called 

grating lobes. To avoid these grating lobes, one chooses 2L> 1, so that the above equation 

d 

will have only one maximum corresponding to n=0. Therefore in this paper, we shall be 
primarily considering values of d up to 1A for most results. 



14 



The nulls of S r (<t>) occur when the numerator alone goes to zero, i.e., whenever 

e jNk c d cos <t> _ ] ? assuming k d cost)) * 27tn. This means cos <b = + , where m is 

Nd 

2 

any integer. As such, for _< 1, there are several solutions for ()> and m, which implies the 

d 

presence of multiple nulls, which might dramatically affect the SINR if they coincide 
with the direction of the desired source. The above example assumed equal weights and 
the absence of noise and mutual coupling. If these were taken into consideration, the extra 
lobes would, in general still be present, albeit shifted. 

D. EXPRESSIONS FOR W opt AND SINR 

With the results of the first three sections, we are now ready to proceed to derive 
the expressions for the optimum weighting vector (w opt ) and the signal-to-interference- 
plus-noise ratio (SENR). 

To begin with, let the desired signal be of the form Ao(t) V 0 OC and the interference 
signal be of the form A k (t) Vk° c where 



A 0 (t) = aoexp[j(cot + p 0 )] 


(23) 


Ak(t) = a k exp(j(cot + p k )] 


(24) 


(V 0 OC ) n = V 0 exp[i(ko d (n-1) cos <f> 0 )] 


(25) 


(V k oc )„ = V 0 exp(j(ko d (n-1) cos fe)] 


(26) 


ao , a k = amplitudes of signals 




p 0 , P k = phases of signals 





15 



Then the input signal to the adaptive processor is 
V,= Z c V oc + N(t) 



= Z c (Ao V 0 0C + X A k V k oc ) + N(t) (27) 

*=i 

1 . Optimum Weight w opt 

We now attempt to find a convenient expression for the covariance matrix <f>. 
Using equation (19) : 

o = <v t v, H > 

= < (Zc V oc + N) (Z c V oc + N) h > 

= Z c < V oc V ocH > Z C H + < N V ocH > Z C H 
+ Zc < V oc N h > + < N N h > 

The following assumptions are made : 

i) The noise voltage at each port is independent of the each other and of the 
signals; its pdf is assumed to be zero mean Gaussian. 

ii) The open circuited desired signal is assumed to be independent of the open 
circuited interfering signals. 

iii) The open circuited interfering signals are assumed to be independent of 
one another. 

iv) The amplitude, phase and angle of arrival of the signal from each source 
vary independently of each other. 

It is shown in Appendix A that this eventually reduces to : 

o = V~ w 2 + Vt Zc Co Zc H + l~; Zc X q k C k Z c H 

= V£ [U + Zc^oCoZc 11 +Z c S^ k C k Z c H ] (28) 
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where = mean square noise voltage at each port 

U is the N x N identity matrix 
Z c = as defined in eqn ( 1 2) 



[C k ]mn= e -> k ° d(ln-n)C0S ^ 



9J k 0 d(m - n)sin <j) k r k ^ 
Rk _ 

k 0 d(m - n)sin <> k r k 
Rl- 



, k = 0,1,2,. ..K (29) 



= input signal (desired source) to thermal noise ratio = 



Vo<a d 2 > 




Igk = input interference (k th interfering source) to thermal noise ratio = — - — — - 

Vn 

Ji (.) = Bessel’s function of the first kind of order one 
It is also shown in Appendix A that : 

S = <A r * V t > 

= V^Z c C s (30) 



where 

2J ( k 0 d(m-l)sin^ 0 r 0 N 

rc 1 =e JM(m-l)cos<D 0 l Rq _ I 

k 0 d(m - l)sin (|) 0 r 0 

R 0 

Substituting eqns (30) and (3 1 ) into eqn (18): 

"opt = d) s 

= ao [U + Z c ^ 0 C o Z C H + Z c X ^ k C k Z C H Y 1 Z c C s 
where a<, is a trivial constant factor and is subsequently dropped. 



(31) 



(32) 
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2 . 



Signal-to-Interference-Plus-Noise Ratio 



With the weight thus calculated, we can proceed to evaluate the SINR, which 
gives a measure of the effectiveness of the array : 



SINR 



Power Received from desired signal 
Noise Power + Interference Power 



Pq 

K 

P N +IPk 



k=l 



(33) 



where, from Appendix B, 

P 0 = < (w H Z c A 0 V o 0£ ) (w H Z c Ao V 0 0C ) H > 

= Vj ^w H Z c C 0 Z c H w (34) 



and 

Pk = < (w H Z c A k V k oc ) (w H Z c Aic V k 0£ ) H > 

= VJ w H Z c ^C k Z c H w (35) 

P N = < (w H N) (w H N) h > = V w 2 w H w (36) 

Therefore, SINR = g 0 w"Z c C o Z c "w (37) 

w H w + w H Z c I^ k C k Z c H w 
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III. RESULTS, OBSERVATIONS AND DISCUSSIONS 



This chapter discusses the results obtained and presents explanations of these 
observations. We first start by asking whether mutual coupling is at all significant and 
then proceed to look at how separation distance, noise level, scatterer radius and number 
of interference sources affect the performance of the array. 

A. SPATIAL CROSS CORRELATION 

The first question we need to ask is whether the effect of mutual coupling is at all 
significant to an adaptive array? The importance of mutual coupling to diversity arrays 
has been addressed in [Ref. 7], The parameter of interest is the spatial correlation of the 
field. It is an indication of the correlation of the voltage received between the elements of 
an array. One of the benefits of an array is to provide diversity; this means that while one 
element may be in the deep fade of the received signal pattern, hopefully another would 
still be in a position to receive the signal. For this to be true, the elements have to be 
relatively uncorrelated - hence a small value of spatial correlation means greater 
diversity. 

Consider a two element array, distance d apart. In the absence of thermal noise, 
the voltage equations from (14) becomes : 

V, = Zcn V, 0C + Zc, 2 V 2 0C 
V 2 = Zc 2 1 V i oc + Zc 22 V 2 ° c 

Then the spatial cross correlation of the received electric field is defined as : 

p = < V, V 2 * > = Zc ii Zc i 2 * <Vi 0C V i° c *> + Zc ii Zc 22 * <Vi 0C V 2 0C *> 

+ Zc 12 Zc 2 r <V 2 0C V, oc ’> + Zc 1 2 Zc 22 * <V 2 0C V 2 0C ’> (38) 
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<V,“ V, <v > = <V 2 ot V 2 ‘* *> = 1 




(39) 


and 


<vr vr > = <V 2 “ V, <x -' *> = J„ (k„d) 


|Ref. 4] 


(40) 


where 


Zcil=Zc22 , Zci 2 = Zc 2 l 




(41) 



J 0 (.) = Bessel’s function of first kind, order zero. 



Thus 

p = (ZcnZc ,2 + Zc 12 Zc ii ) + (Zc u Zc 11 + Zc 12 Zc 12 ) Jo(k„d) (42) 
Extracting the common factor and substituting for Z c , the spatial cross correlation 
of the electric field for a two element array terminated in Zl. = Zn* is given by [Ref. 4] : 



P — pot 



R 



12 



R, 



1 + 



|Z|2 | 

AR 2 n 



\- \ 



(43) 



v j 

where p,*- = J c (kod) is the correlation in the absence of mutual coupling, and R n and R, 2 
are the real parts of Zn and Z 12 respectively. Fig 3.1 below shows the spatial cross 
correlation vs separation distance d between two half-wavelength dipoles. 



Spatial Cross Correlation of a 2-element array 




Fig 3. 1 Spatial Cross Correlation of a 2 element array 
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As expected, at longer separation, the effect of mutual coupling is diminished and 
the two curves converge. However, there is a marked difference for d<X, which indicates 
the significance of including mutual coupling. 

Another way to interpret spatial cross correlation is to examine the normalized 
spatial cross correlation (Fig 3.2). For reference, we have indicated the correlation values 
of 0.7 and 0.3 (Table 3.1). It is seen that for this example the effect of mutual coupling is 
to decorrelate the signals initially but correlate them more for larger values of spacing d. 

Normalised Spatial Cross Correlation of a 2 elemert array 




Fig 3.2 Normalized Spatial Cross Correlation Electric Field in a 2 element array 





Values of d/a. 


Without mutual coupling 


With mutual coupling 


p < 0.3 


> 1.1201 


> 3.0701 


p < 0.7 


>0.1801 


> 0.0901 



Table 3. 1 Values of spacing d for p <0.3 and 0.7 
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While both curves fall rapidly below the upper level of 0.7, there is a large 
difference in the lower level of 0.3; without mutual impedance, the arrays have minimal 
spatial cross correlation (<0.3) for d>1.12X, whereas with mutual coupling, this occurs 
only after d>3.07L What this means is that in the presence of mutual coupling, the array 
diversity will be generally reduced. 



B. EFFECT OF MUTUAL COUPLING ON WEIGHTING VECTOR 

Another effect of mutual coupling is in the resultant weights. Using (18), we 
observe that we get different weighting functions when considering and ignoring the 
effects of mutual impedance (Table 3.2). The weights are normalized such that w H w = 1. 



Number of elements 
Separation between elements 
Number of interference sources 
Direction of desired source 
Directions of interference sources 

Scatter radii 


N = 4 
d = A/4 
K = 3 
£ = 90° 

<p k = 30°, 60°, 120° 
r 0 = r k = 0 


Consider mutual coupling 


Ignore mutual coupling 


Deterministic case 


w (normalized) 


Wo(normaHzed) 


w d (normalized) 


0.2624-0. 1203i 


0.391 l+0.0261i 


0.3300 


-0.4627-0. 2327i 


-0.3631 -0.463 li 


-0.5357-0.3227i 


0.4330+0.5600i 


0.2335+0.5445i 


0.427 5+0. 4564i 


0.2102-0.321 li 


0.1704-0.3530i 


-0.0689-0.3227i 



Table 3.2 Weights for adaptive MMSE and deterministic cases 



Table 3.2 also shows the weights Wd calculated in the deterministic case, that is, 
where it is assumed that the directions of the interference sources are known and the 
weights calculated such that nulls are steered in the directions of the interferences. Since 
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the weighting veetors are all different, we would expect different beam patterns using 
these 3 solutions, as shown in Fig 3.3 below. Values of the corresponding SINR arc 
indicated in the inset. 

Note : In Fig 3.3 and subsequent plots, the blue curve corresponds to adaptive (MMSE) 
solution using cqn (18) , with mutual coupling considered. This means that the off 
diagonal terms of Z are non-zero and are given by (6). The red curve corresponds to 
adaptive solution also given by (18), this time ignoring mutual coupling (Z in (13) is 
diagonal). The green curve corresponds to the solution in the deterministic case, 
calculated without mutual coupling. Moreover, the direction of the desired source is 
maintained at <p 0 =90° (unless otherwise specified) and is denoted by a blank circle on the 
abcissa. Pink asterisks (*) denote the directions of the interference sources. 

Weights normalised N= 4, d= 0.25 a, K = 3, rjf\ = 0, r k /R = 0, = 20dB, 2 k = 20dB 




Angle of Arrival (In degrees) 

Fig 3.3 Beam pattern for adaptive solutions considering and 
ignoring mutual impedance and deterministic weights. 

All weights normalized such that w 11 w=l, beam pattern not normalized. 
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In Fig 3.3, it is seen that although nulls are steered in the directions of the 
interference sources in the deterministic ease (green), it does not provide the best SINR. 
This is because in the other cases (adaptive solutions - blue and red), the gain in the 
direction of the desired signal is higher than the deterministic case (which does not take 
into account direction of desired signal), thereby contributing to a higher desired signal 
power received, P Q in (34). This thus gives rise to higher SINR. 

Fig 3.3 shows the beam pattern (un- normalized) using weights calculated. 
However, it is often the normalized beam pattern that is more suitable for comparison 
between the three solutions; in Fig 3.4 and subsequent graphs, the beam patterns are 
normalized to the maximum gain. 

N= 4, d= 0.25A., K = 3, r Q /R = 0, r k /R = 0. £ Q = 20dB. = 20dB 




Fig 3.4 Beam pattern (normalized to max gam) lor adaptive solutions considering and 
ignoring mutual impedance and using deterministic weights 



Fig 3.3 and Fig 3.4 are essentially the same beam patterns, though in Fig 3.4 the 
adaptive, no-mutual-coupling solution (red) seems to have a higher gain in the direction of 
the desired source {(p 0 = 90°) than the mutual coupling case (blue). As such, one might be 
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lulled into expecting the red curve to have a higher S1NR; it should be borne in mind that 
the graph shows relative gains only. For a true value of the SINR values, one should 
always refer to the values displayed in the inset rather than to the heights of the curves. 

The previous two sections show that mutual coupling does have an effect on the 
performance of adaptive array antennas. This effect obviously diminishes as the physical 
separation between the elements increases. 

C. EFFECT OF ELEMENT SEPARATION 

An interesting comparison would be how the SINR is affected by varying the 
separation distance d between elements. Fig 3.5 below shows the SINR vs d plot for a 4- 
element array antenna in a 3 interference situation (no scatterers). 



N= 4, K = 3, r o /R = 0. r k /R = 0, = 20dB, = 20dB. (j) k = (30°, 60°, 1 20°) 




Fig 3.5 SINR vs d 
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It is seen that the SINRs with and without mutual coupling do show differences for 
small values of d (less than IX), sometimes by as much as 3dB (e.g. at d=().75X). At 
higher separations, they are about the same, as expected since mutual coupling effect is 
diminished. In all cases, the SINR using the adaptive solution (no mutual coupling) is 
consistently higher than that using the deterministic weights. At the same time, there are 
differences between the adaptive solutions with and without mutual coupling (blue and red 
curves in Fig 3.5). The exact difference is difficult to quantify as it is case specific: in Fig 
3.5 above, the situation is for a 4 element array in a 3 interference environment {<j) k = 30°, 
60° and 120°) with = 2()dB. 

An interesting feature is the fact that the SINR drops drastically near d=1.15X. 
This is shown in Fig 3.6 below. 

N= 4 . d= 1 .1 5X.. K = 3 , r o /R = 0, r^/R = 0 , = 20dB. ^ k = 20dB 




A closer examination reveals that this drop is due to the rather low angular 
response of the array in the direction of the main source. Therefore, in practical design of 
array antennas, it would be wise not to exceed d=X/2 unless we impose an additional 
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constraint that the pattern not go to zero in the direction of the desired source. In this 



paper, we shall continue to study the cases for d up to IX. 



D. EFFECT OF NOISE 

From (18) and (14), we have : 

Wopt = O' 1 S 

= [ U + Zc Co Zc 11 + Z. 2 ^ C k Zc 11 r 1 Zc C s 
and V, = Z c V 0C + N(t) 

In the situation where £,<>» 1, » 1, we get IV 0 *! » IN(t)l so 

W„p, « [ Ze & C„) Ze H +Ze(i:^C k )Z c H ]- 1 Z c C s (44) 

and V.^ZcV 00 (45) 

The resultant signal is 

S r = Wop, H V, 

» {[Zc^CoZc 11 +Z £ I^C k Z c H r 1 Z c C s } H Z c V oc 
= C s H Z c H {[Z c ^oC 0 + Z^ k C k 
= C S H Zc H {(Zc 11 )- 1 & C 0 + I Ct)' 1 Zc' 1 } : H Zc v~ 

= C S H Zc 11 (Zc 1 ) H [&> Co + Z 5k Ck)' 1 ]® [(Zc 0 )' 1 ] 0 Zc V" 
= C 0 (Zc-'Zc) 0 [(£> Co + Z ^ Ck)' 1 ] 0 [Zc H (Zc 0 )- 1 ] 0 V~ 

= C S H [(£, C 0 + 1Z* Ck)- 1 ] 0 V 00 (46) 
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which turns out to be independent of Zc. Since the effect of mutual interference is in the 
impedance matrix Zc, it is predicted that when E, 0 , £ k is large (>40dB), the effect of mutual 
impedance is reduced and the two adaptive solutions (blue and red in Fig 3.7) converge. 
(Recall that £, 0 = input signal (desired source) to thermal noise ratio and = input 
interference (k Ul interfering source) to thermal noise ratio ) 

N= 4. d= 025X. K = 3. r o /R = 0. r^R = 0. = 20dB, = 20dB 



(a) 





Fig 3.7 Comparison of the effect of noise 



Fig 3.7a shows the beam patterns using the different weights when E, (> = £ k = 2()dB. 
However, when signal-to-noise-ratio £, 0 and interfcrencc-to-noise-ratio £, k are increased (to 
40dB) as in Fig 3.7b, one can easily observe that the solutions (with and without mutual 
coupling) converge, as predicted. It is important to note that though the normalized beam 
patterns in Fig 3.7b are similar, the weights and terminal voltages with and without mutual 
coupling are not the same. This is evident in the fact that the SINR values with and 
without mutual coupling differ by as much as 4dB in Fig 3.7b. This is because the noise 
levels in the three cases arc different relative to the output signal. 
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E. EFFECT OF SCATTERER RADII r„ , r k 

Recall that in our model outlined in Section 2A that the desired and interference 
sources are at distance R k from the array. In addition, each source is embedded in a 
scattering environment, in which the scatterers are assumed to be distributed uniformly 
around each source over a circular region of radius r k (k=0,l,....K). Here we take R k = 1 



km and vary the values of r k (k=(),l,....K); ultimately, it is the ratio b k =— that is 

R 



significant in the calculations. 

It was found that the effect of the scatterer radii r Q and r k depend on the input 
signal to noise ratios and ^ k . First we compare the effect of increasing scatterer radius 
r„ while maintaining r k = 0. For illustration, a 4-element array with A/2 spacing is used in a 
3 interference situation (<p k = 30°, 45°, 60°). Fig 3.8 below shows the resultant normalized 
beam pattern. 



(a) 



(b) 



N= 4, d= 0.5X, K = 3, r Q /R = 0. r k /R = 0, = 10dB, ^ k = 10dB 




N= 4, d= 0.5X, K = 3, r o /R = 0.1 , r R /R = 0. % Q = 1 0dB, = 1 0dB 
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In Fig 3.8a, there are no scatterers around the desired source (r c = 0) whereas in 
Fig3.8b, scatterers are present (ro/R = 0.1). We can see that the weights and (hence) beam 
patterns are changed slightly, and that with scatterers (Fig 3.8b) the resultant SINR values 
are slightly reduced. The light blue column shows the extent of the angular spread of the 
desired signals; hence the area in the column under each curve would represent the 
amount of power arriving from the source. In the calculation of SINR, this area is divided 
by the angular spread to give an average value of the signal. Whether this average value is 
greater or smaller than the SINR in the no-scatterer case obviously depends on the shape 
of the curve. In Fig 3.8 above, it is observed that though the average SINR values for 
adaptive solutions are smaller, the deterministic solution actually shows an improvement in 
the averaged case. 

Nonetheless, it is important to note that the adaptive solutions (blue and red) still 
give greater SINR than the deterministic case (green). This is in accordance with our 
conjecture that the adaptive MMSE solution also leads to the Maximum SINR solution. It 
has been shown in [Ref. 4] that the SINR without mutual coupling and no scatterers 
around the desired source is always better than the deterministic case. No such conclusion 
can be drawn when mutual coupling is taken into account. 

Next we investigate the situation when and are increased to 20dB (Fig 3.9); 
comparing the 4 cases (in Fig 3.8 and Fig 3.9), the deterministic case remains the same 
since the deterministic weights do not depend on r 0 , r k , nor % k . The adaptive solutions 
change accordingly - most surprising is the fact that for r</R = 0.1, r k /R = 0, ^ = 

20dB (Fig 3.9b), the adaptive solutions (both with and without mutual coupling) give 
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SINR values smaller than the deterministic case. This is in contrast to the no-seatterer ease 



where the adaptive (MMSE) solutions always give the maximum SINR. 

N= 4, d= 0.5X, K = 3. r /R = 0, r k /R = 0, £, = 20dB, ^ = 20dB 



(a) 



(b) 





Next we examine the effect of increasing the seatterer radius r k while maintaining 
r« =(). In Fig 3.10 below, the columns in magenta represent the angular spread of the 
interference signal from the interference sources. The increment of r k again reduces the 
SINR but not as much as in the case of scatterers around the desired source. The decrease 
in SINR is expected since scatterers increase the directions from which the interference 
signals can now reach the array, thereby increasing the interference power received. In this 
case the increment of r k does not change the fact that adaptive solutions still give greater 
SINR than the deterministic case. 
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(a) 



(b) 



N= 4, d= 0.5/., K = 3, r o /R = 0, r R /R = 0, = 20dB, £ k = 20dB 





Fig 3. 10 Effect of scatterer radius r k (<;„ and % k = 2()dB) 



F. EFFECT OF NUMBER OF INTERFERENCE SOURCES K 



Another variable parameter is K, the number of interference sources. In Fig 3. 1 1 
below, a 4-element array (N=4) is placed in an environment where there are no scatterers 
around any of the sources (r„ = r k = 0). 

N= 4, d= 0.5/., K = 3, r/R =0. r k /R = 0,\ o = 20dB, § k = 20dB 



(a) 





Fig 3.1 1 Effect of Number of Interference Sources K (K>N-1) 
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In Fig 3.11a, 3 interference sources (K=3) have been specified in the directions 
<p k = 3()°, 60°, 120°. The adaptive solutions and the deterministic solution are similar and all 

steer nulls in the directions 30°, 60°, 120° that effectively cancel out the interference 
sources. The resultant SINRs for all 3 cases are therefore very close. In Fig 3.11b, an 
additional interference source has been added (at <p k = 140°). With the deterministic 

weights, a maximum of N- 1 nulls can be chosen; this means we can (arbitrarily) choose 
only 3 of the 4 interference sources to cancel. With the adaptive solutions specified by 
(18), the weights are altered to give a different beam pattern as shown. Though none of 
the nulls is exactly in the direction of the interference sources, the ‘spreading’ of the 
pattern is such that the resultant SINR is greater than that given by the deterministic case. 

In Fig3.12 below, the same situations are compared, this time with scatterers 
around all sources. 



N= 4, d= 0.5X, K = 3, r o /R = 0.1 , r k /R = 0.1 , % o = 20dB, = 20dB 





Fig 3. 12 Effect of Number of Interference Sources K (K>N-1) 
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In Fig3.12a, wc note that the deterministic SINR is better than for the adaptive 
eases (as discussed previously in Section E). However, with the addition of a fourth 
interferer at <p k = 140°, the adaptive cases actually perform better than the deterministic 

solution in terms of SINR. This emphasizes the llexibility of the adaptive solution in 
encountering more than N-l intcference sources. 

At the other extreme, we would like to examine the situations where the number of 
inteferers K is less than the number of array elements N. To illustrate, we consider N=6 
and 8 in an environment with only 3 interference sources. Fig 3.13 below shows the 
resultant normalized beam patterns. 



N= 6, d= 0.5X, K = 3, r o /R = 0. 1 , r k /R = 0.1 , = 20dB, 4 k = 20dB 





Fig 3.13 Effect of Number of Interference Sources K (K<N-1) 



It is seen that when K < N-l (“there are more than enough array elements to go 
around”), the SINR using the adaptive solutions are improved. Table 3.3 below shows 
these values for N=4,6,8,10 and 30 (K=3 in all cases). 
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K = 3 d = X/2 Z^iL-oi = 20dB 

R R 


N=4 


N=6 


N=8 


N=10 


N=30 


SINR (with mutual coupling) 


15.4dB 


15.9dB 


17.8dB 


18.2dB 


18.5dB 


SINR (ignore mutual coupling) 


l5.2dB 


15.7dB 


17.4dB 


17.7dB 


18.0dB 


SINR (deterministic case based on N=4) 


16.9dB 


16.9dB 


16.9dB 


16.9dB 


16.9dB 



Table 3.3 Comparison of SINR 



As is obvious from the table, as the number of array elements (N) increases, the 
SINR increases as well for the adaptive cases. Of course, the deterministic case does not 
change since it is dependent only on the direction of the interferers. Initially, we note that 
SINR (deterministic) is greater than SINR(adaptive) - this was discussed in Section 3E as 
being due to scatterers around the desired source. However, for higher values of N (N> 8 
in this example), we find that the SINR of the adaptive solutions actually increases to 
greater than that offered by the deterministic case. Beyond a certain value of N, this 
increment becomes insignificant or non-efficient (c.f. SINR for N=10 and 30). 
Nonetheless, this is another indication of the ability of the adaptive algorithm in making 
full use of its elements to find the best solution. 

G. ANGULAR RESOLUTION OF SOURCES 

Our next area of interest is in examining if the use of an adaptive array is effective 
when the various sources are closely spaced. Consider the situation where one of the 
interferers is very close to the desired source : (j) 0 = 90°, tp k = 90°-2.85° = 87.15° (2.85° 
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V V 

was chosen as it is the rms value of angular spread for — =— = 0. 1 as in eqn (5) ). This 

R R 



is shown in Fig 3.14. 



N= 4. d= 0,5/.. K = 3, r /R = 0, r k /R = 0,1 = 10dB. ^ k = 10dB 



(a) 



(b) 





3.14 Effect of Interference Source close to Desired Source (^ D and ^ k = lOdB) 



Fig 3.14a shows the situation where there are no scatlerers around any of the 
sources. It is seen that the close proximity of a null near the desired source lowers the 
SINR in the deterministic case drastically, whereas for the adaptive solutions, the nulls are 
actually steered away so as to give a significantly better SINR. In Fig 3.14b, the presence 
of scatterers does affect the performance of the adaptive solutions, but it is seen that the 
adaptive solutions still give better SINR values. 

Fig3.15a and b below shows the same situations with and £, k increased to 40dB. 
As discussed in Section D, high values of £, 0 and mean that the solutions all converge to 
the deterministic case. However, in the presence of interference close to the desired 
source, the solutions are different. More significantly, the SINR from the adaptive 
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solutions in this case is lower than that from the deterministic case - again, this is due to 



the presence of scattcrcrs around the desired source. 



N= 4, d= 0.5A, K = 3, r /R = 0, r. /R = 0, £ = 40dB, l. = 40dB 
’ ’ ’ o ’ k ' r ^k 





3.15 Effect of Interference Source close to Desired Source (£,„ and £, k = 40dB) 



What happens when an interference source, initially remote from the desired 
source, is moved close to it? In Fig 3.16a below, <j) k = 30°, 45° and 60° while in Fig 3. 16b, 

the interference source at <p k = 45° is moved to 87.15°. 
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N= 4, d= 0.5*., K = 3, r /R = 0.1 , r./R = 0.1 , £ = 20dB , L = 20dB 



(a) 



(b) 





60 80 100 120 
Angle of Arrival (in degrees) 



160 180 



Fig 3. 16 Effect of Interference Source close to Desired Source (E,„ and = 20dB) 



In a previous section, it was shown that when there are scatterers around the 
desired source, the deterministic solution gives a greater SINR than the adaptive solution, 
as shown in Fig 3. 1 6a. However, this is reversed when one of the interference sources is 
moved close to the desired (Fig 3.16b). This shows that when an interference source is 
close to the desired source, the adaptive solution performs better than the deterministic 
one, even when there are scatterers around the desired source. 
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IV SUMMARY AND CONCLUSIONS 



A. SUMMARY 

Through the senes of results obtained, the effect of mutual coupling on the 
performance of adaptive array antennas was studied by varying several parameters. These 
can be summarized as follows: 

1. Effect of Mutual Coupling on Spatial Cross Correlation 

The effect of mutual coupling is to reduce the spatial cross correlation initially 
(for d < IX) but to increase it for larger values of inter-element spacing d. 

2. Effect of Separation Distance d 

For d>l/ ) i, there is practically no difference on the SENR obtained with and 
without considering mutual coupling. At lesser separation, however, the difference is 
significant and can be as much as 3dB. 

3. Effect of Noise 

When the signal-to-noise-ratios £, 0 and ^ are high (>40dB), it was found that the 
solutions obtained with and without considering mutual coupling are very close, and 
converge to the deterministic solution (provided number of interference K < number of 
array elements N). 

4. Effect of Scatterer Radii r 0 and r k 

When the scatterer radii r Q and r k are increased, in general, the SINR of the 
solutions are reduced, This is expected as it represents uncertainty in the directions of 
arrival of the sources. 
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When the scatterer radius r 0 is increased, at a critical £#, ^ level, the adaptive 
MMSE solutions cease to give Maximum SINR. This is in contradiction to the conjecture 
that MMSE and Maximum SINR solutions are equivalent [Ref. 6: p4 1 ]. 

On the other hand, increasing the scatterer radius rk alone does not ‘invert’ this 
trend of MMSE solution giving max SINR. 

5. Effect of Number of Interference Sources K 

When the number of interferers K > the number of array elements N, the 
deterministic solution can do no better than pointing N-l nulls towards the interferers. 
The adaptive solutions, however, seeks to minimize the error and generally gives a 
greater SINR than the deterministic case. 

When the number of interferers K < the number of array elements N, it was found 
that the adaptive algorithm makes use of this fact to improve the SINR. Hence even with 
scatterers around the desired source (as in the previous section), with a sufficient number 
of array elements, the adaptive solution will give an SINR greater than that in the 
deterministic case. 

6. Angular Resolution of Sources 

In cases where one of the interferers is very close to the desired source, it was 
found that adaptive solutions perform better than deterministic one. This includes the 
situation where there are scatterers around the desired source and demonstrates the 
excellent angular resolution offered by adaptive arrays. 
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B. CONCLUSIONS 



We have examined the effect of mutual coupling on the performance of an 
adaptive array in a realistic communications environment, where the desired and 
interference sources are embedded in an area of scattering objects. Expressions for the 
MMSE solutions and SINR were derived for this scenario. These were subsequently 
used to show the effect of mutual coupling, in particular in array where the element 
separation was less than 1^.. The effects of the other parameters in the system were also 
discussed, including conditions under which the MMSE solution does not give rise to 
Maximum SINR. 
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APPENDIX A : EXPRESSIONS FOR O AND S 

A. EXPRESSION FOR O 
From (15), 

covariance matrix <E> = <V t V t H > 

= < (Z c V oc + N) (Z c V oc + N) h > 

= Z c < V oc V oc H > Z C H + Z c < V oc N h > (Al) 

+ <NV ocH >Z c H + <NN H > 

Now, consider the terms separately : 

(i) Second and third terms : 

Assuming noise is zero mean and uncorrelated with signal, 

then < ao e _j(co ° t+|3 ° } N(t) > oc < N(t) > = 0 

(A2) 

so < V, oc N H > = < V, oc > <N h > = 0 = < N V t ocH > (A3) 

(ii) Fourth term : 

Assuming noises at the ports are uncorrelated with each other 

< N N h > mn = 8 mn v£ 

(A4) 

where V N = root mean square noise voltage 

(iii) First term : 

< V, oc V t 0£ H > = < [ao e j(<ot+po) V 0 0C + Z a k e i( ““ +pk) V k oc ] 

. [ao e' j(CM+p0) V 0 0C H + Z a, e ' j((0t+pi) V, oc H ] > 

= < ao 2 V„ oc (V 0 0C ) H > + Z < ao a k e j(pk ' Po) V k oc (V 0 ° c ) H > ( A5) 
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+S < a,, a k d (p0 ' pk) V 0 0C (V k oc ) H > +IZ < a k a, e i((3k ' pl) V k oc (V, oc ) H 



> 



Now <a 0 2 V 0 0C (V 0 0C ) H > = <a 0 2 ><V 0 0C (V 0 0C ) H > (A6) 

(assuming magnitudes of narrowband signals are uncorrelated to the angles of arrivals) 

Also, assuming the interfering signals are uncorrelated with the desired signal, 
then <aoa k > = <ao> <a k > = 0 
since the phase (5 0 is assumed uniform over (0, 2k). So 

< ao a k e i((3k ' p0) V k oc (V 0 0C ) H > = <ao> <a k > <e iPk > < e i(5o > <V k oc xV 0 oc H > (A7) 

= 0 

and < a k a, e i(pk ' pi) V k oc (V, oc ) H > = 5 kl <a k 2 V k oc (V k oc ) H > 

= 5 kl <a k 2 xV k oc (V k oc ) H > (A8) 



Now consider <V k 0C (V k oc ) H > = V 0 2 < e jk 0 d(m-n)cos<!> k > 

= V 0 2 < e jk 0 d(m-n){cos cos a k -sin <ji k sin a k } > 

.k 0 d(m-n) - / 2 Z2 ■ t 

J {cos<t> k yb k -x -sin <j) k x} 

= V 0 2 < e bk > 

where x = —sin a k (A9) 




j k°d(m 0 ) j CQs ^ ^b k 2 -x 2 -sin <j) k x } 

bt. 




A 

, fl k^dCm-n) — J “ 
=V 0 “ j ^ cos( — — -xsin (ji^le 



•k 0 d(m-n) - 






-cos 






-X 2 dX 



2 -x 2 



— Vl-x 2 dx 

K 
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Let us now assume that b k » VlO =>b k 2 » 10, so that-^/b k ~ -x 2 = b k , then 

<v k 0c (V k 0C ) H > * V 2 e jk 0 d(m-n)cos^ k 2 f 1 Vl - X 2 cos(- X sin <j> k )dx 

7t J ~l Ri- 



2J i 



_ y 2 e jk 0 d(m-n)cos4> k 



- r k 

(k 0 d(m-n)stn 4> k — ) 
K k 



• — 

(k 0 d(m-n)sin<t> k — ) 
R k . 



= V 0 2 e jM( n '-»)cosfc fm _n (kodsin4)kirkiRk) 



(A10) 



where f m _ n (k 0 d sin <|> k ,r k , R k ) = 



2J r k 0 d(m- n)sin([> k r k ^ 

Rk _ 

k 0 d(m-n)sin<)) k r k 
Ri- 



(All) 



and Ji = Bessel’s function of the 1 st kind, of order 1. 



Substituting (A9) into (A4), 



< V t oc V t ocH >mn = V 0 2 { <a 0 2 > e jk o d(m - n)cos<t, °f m _ n (k 0 dsin(t) 0 ,r 0 ,R 0 ) 



+ Z<a k 2 > e jk o d ( m-n ) cos ^ f m _ n (k 0 d sin (Jtk , r k , R k ) } 



Therefore, < V t oc V t oc H > = V 0 2 { <ao 2 > C 0 + 2 (<a k 2 > C k ) } (A 1 2) 
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where [CJ™ = e jk “ d(m ' n)coi 



k 0 d(m-n)sin<|> k r k 

2h R 

^jc 

k 0 d(m - n)sin <j) k r k 

Ri- 



(A13) 



Finally, O = Z c < V, oc V, oc H >Z C H + < N V, oc H > Z C H + Z c < V t oc 
= V 0 2 + V 0 2 Z c { <ao 2 > C„ + I (<a k 2 > C k ) } Z C H 

= vg + V^ Z c & Co) Z C H + v£ Z c (X q, CO Z C H 
= V§ [U + Zc(^oCo)Zc H + Zc(X^C k )Z c H ] 



where ^,= 





> 




is the input signal to thermal noise ratio 



N H > + < N N H > 



(A 14) 
(A 15) 



* V 0 2 <a k 2 > 

Sk = 



V 



N 



is the input interference to thermal noise ratio 



(A 16) 



B. EXPRESSION FOR S 

From (16) S = <A r H V t > 

where A r (reference signal) takes the form ao V 0 e^ to ° t+ ^ 0 ^ (A 17) 

then S = < ao V 0 e _j(0) ° t+Po) {Z c [ao e j(c0 ° t+Po) V 0 0C + J a k e j(o) ° t+Pk } V k oc ] + N(t) } > 

1 

= V 0 Z c < ao 2 (V 0 0C ) + X (ao a k e*®* ' Po) V k oc ) > (A 1 8) 

+ V o <aoe _j((0 ° t+Po) N(t)> 



<V 0 0C > n 



= V < e jk 0 d(n-l)cos<t) 0> 
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- v 0 e ^ 0<1 ( n - | > cos< t' 0 f n- | ( k 0 dsi n <t> 0 .r 0 . R ) 

Lastly, < ao ak e^ (pk ' p0) Vk° c > = < a« ak > <e i(Pk ‘ Po) > <V k 0C > = 0 
where < ao ak > = 0. 

Therefore, using (A20), (A23) and (A24), 

S =V 0 2 <a 0 2 >Z c [f n _ 1 (k 0 d,r 0 ,R 0 )] 

= ^ZcC s 

where 



2J t (k 0 d(m - 1) sin <() 0 ~~) 

[C s ] m =e jk ° d(,n_1)cos< *’ 0 K ° 

(k 0 d(m-l)sin <j) 0 



(A 19) 
(A20) 



(A21) 



(A22) 
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APPENDIX B : EXPRESSIONS FOR P 0 AND P 

From (32), we have 

P 0 =<(w h Z c A„V o oc )(w h Z c A o V o 0c ) H > 

= <0o 2 > w H Z c <Vo OC Vo OC H > Zc H W 
= w H Z c Co Z C H w 



From (33), we have 

P k = < (w H Z c A k V k oc ) (w H Z c A k V k oc ) H > 
= <a k 2 > w H Z c <V k 0C V k oc H > Z c H w 

= v£ w H Z c (Sg k C k ) Z c H w 



where V N = root mean square noise voltage 
, ^ k as defined in (A 16) and (A 17) 
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